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ABSTRACT 


An analytical solution is developed for the determination of the 
stresses and displacements in a unidirectional fiber -reinforced composite 
containing an arbitrary number of broken fiber as well as longitudinal 
yielding and splitting of the matrix. 

The solution is developed using a "materials -mode ling" approach which 
is based on a shear-lag stress transfer mechanism. The equilibrium equation 
in the axial direction gives a pair of integral equations which are solved 
numerically. 

Excellent agreement is shown to exist between the solution and 
experimental results for notched unidirectional boron/aluminum laminates 
without splitting. For brittle matrix composites (ie. epoxy) equally good 
results are indicated for both matrix yielding and splitting. 

For yielding without splitting the fracture strength is found to 
depend on crack length while for large splitting it is crack length 
independent. 


/l/f < d ~ /zTd 6 <, 
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INTRODUCTION 

An attempt is made in this study to develop an analytical model capable 
of predicting the characteristic strength and fracture properties of a 
unidirectional composite laminate . The investigation considers a two- 
dimensional region containing an arbitrary number of broken fibers as well 
as longitudinal matrix damage as shown in Figure 1. The fiber breaks lie 
along a transverse line and therefore represent a notch. Damage to the 
matrix originates at the ends of the notch, i.e., in the region between 
the last broken and first unbroken fiber and consists of both yielding 
and splitting. Symmetry is assumed as indicated and only the first 
quadrant of the region is shown in Figure 1. 

The matrix is taken to be an elastic -perfectly plastic material and 
the fibers are linearly elastic. Load is transferred from adjacent 
fibers through the matrix by a simple shear-lag mechanism with the shear 
stresses being independent of transverse displacements. The axial fiber 
stress is also independent of transverse displacements and the equilibrium 
equation in the fiber direction reduces to an equation in the longitudinal 
displacement alone, as is typical of shear-lag solutions. 

Similar investigations have been presented by Hedgepeth and Van Dyke 
in [l] and [2] in which only one broken fiber was considered with yielding 
alone in [l] and splitting alone in [2]. The extension to more than one 
broken fiber however is not developed conveniently by the influence function 
technique as suggested in [l] because the broken fiber adjacent to the 
damaged region is not typical of any of the remaining broken fibers . 
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For a given number of broken fibers and. a known applied axial stress, 
taken to be uniform at points remote from the damage, the stress concen- 
tration in unbroken fibers, the extent of the matrix damage and the fiber 
displacements are desired. 

An interesting study is presented by Peters in [3] concerning the 
fracture strength of unidirectional composites which exhibit large matrix 
splitting (Boron/epoxy, graphite /epoxy ) and those such as boron /aluminum 
which have large plastic yielding but little splitting. For the first class 
the fracture strength is independent of crack length while in the second 
crack length dependence is found. This behavior is considered in detail 
using the present model and the extremes of large splitting and large 
yielding with no splitting are predicted accurately. The. model apparently 
does hot have the capability of accounting for those composites which, 
exhibit matrix yielding and small but stable longitudinal splitting. Some 
consideration is given to the reasons for such difficulties and possible 
modifications to give a. more complete model, are discussed. 



3 

FORMULATION 

The laminate is modeled as a two-dimensional region, shown in Figure 1, 
having a single row of parallel, identical, equally spaced fibers with the 
broken fibers being symmetric about the center line and the matrix damage 
occuring between the last broken and the first unbroken fiber. It is 
assumed that the fibers have a much higher elastic modulus in the axial 
direction than the matrix and therefore the fibers are taken as supporting 

°** ax * a ^ stress in the laminate. The matrix supports transverse 
normal stresses and shear stresses. 

Admittedly, most unidirectional composites consist of more than one 
lamina with all fibers in each lamina surely not perfectly aligned either 
through the thickness or, with-in each layer. These variations can have 
a considerable influence on the stress state. For example, in [1] and 
[5] it is shown that the shear stress becomes large as the fiber spacing 
decreases, i.e. 0(l//d) for rigid, fibers where d is the. minimum distance 
between fibers. Local failure may well occur at critical points through 
the thickness in advance of laminate splitting which could give an 
apparent shear stiffhess considerably different from that for the matrix' 
alone. It is assumed that such variations can be accounted for by an 
appropriate choice of the shear modulus G^ and the transfer distance, h. 

It is with this in mind that the following development will be concerned 
with an equivalent lamina where G^ and h are to be determined 
experimentally for any particular laminate. 
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A free -body diagram for a typical element is given in Figure 2, with the 
special condition for the last broken fiber, denoted by n = N, and for 


y £ L that 





where 


<y4>*l , y > & , 

<y - i> s o • .» y < i , and 


( 2 ) 


L equals the total damaged length., % split length, and T q the matrix 
yield stress. 

The equilibrium equations in the longitudinal and transverse directions 


respectively, for all fibers n, with the exception of N and N+l when 
y <_ L, are 


A_ d a„i 

— — ■■■ ■!■■ - ■ + t I . - 1 1 = 0 , and 

t dy 'n+l n 


(3) 


- D J *° 

n+l n 


(M 


For fiber N, y < L, x i N-+1 = " T o <y “ i> > shd the equilibrium equations 


are 
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A- d0 L 

r ■: t <y ~ l> - r = 0 , and 

t dy o 'N 


( 5 ) 


a | - o | + ~ {-t <y - i> + t|J = 0 . (6) 

m N+l ^ H 2 dy ° N 


For fiber N+l, y <_ L, t| n+1 - -T Q <y - &>, and the equilibrium equations 


are 


Ay da. 


___ f 'N+l 
t dy 


+ T N+2 + T o <y “ A> = 0 » and 


( 7 ) 


(8) 


Further simplifying assumptions are now made regarding the stress- 
displacement relations which reduce the number of unknowns from three 
stresses to the two displacements, u^ and v^. Let 


dv 

i n 

a fl F dy - ’ 

n 


( 9 ) 


a | = 2.7 u. - u }/h, and 

m' M n+l n 

n 


(10) 


t| .. = G {v - v }/h . 
'n+l M n+l n 


( 11 ) 


Substituting into the equilibrium equations, the following pairs of 
equations are obtained: 

For all fibers, except N and 'IHi when, y L, 


h V S F 


d v 


n 


G t . 2 
M dy 


{v - 2v + v . } = 0 , and 
n+l n n-1 


( 12 ) 
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E. 


M 


{u — 2u + u . 
h n+l ft n-1 




- v , ) = 0 . (13) 

n+l n-1 


For fiber y <_ L, 


d^v , 

— + v - v„ - 1 <y - i> - 6 , and 

V dy 2 "- 1 K g m 0 


E. 


a. 


£ <W 4 V^’W> * Hr'# ^ 


Vl 1 


(lM 


_ X <y _$,>} = 0 
0 


(15) 


For fiber N+l» y L. 


,2 „ 

d v. 


G,.t ,2 

M dy 


hx 

+ tt * " 7 * + <y “ * 0 y and 

v N+2 N+l G m j 


(16) 


E. 


M r_. o„ + „ V + k 2L {Ji [ v . v 1 

V 2 dy S L N+2 N+1 J 


h {u N+2 “ 2U N+1 


mb x <y ** &> } * 0 » 
O'* 


(IT) 


The shear stress ^displacement form assumed in equation (11) is 

referred to as the shear-lag assumptioh and, as £an be seen above, the 

equilibrium equation in the axial direction is independent of the 

transverse displacement u^. It is then possible to obtain a solution 

for the axial displacement v^, and therefore the fiber stress and shear 

stress , independently of u . Onoe v . is known, the transverse 

n n 

displacement and matrix stress may be obtained from the remaining 
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equilibrium equation. References [5 1 and [ 6 ] consider similar three- and 
two-dimensional solutions, without .matrix damage however, in which the 
shear stress is assumed to depend on the transverse as well as the axial 
displacement and the equilibrium equations do not uncouple as for the 


shear -lag assumption. 

It is the intent of this study to investigate behavior due to broken 
fibers and matrix damage in which the failure criterion for the matrix 
is due to shear alone and the matrix is assumed to be elastic-perfectly 
plastic. In this case, the matrix transverse normal stress plays no role 
and the remaining discussion will focus on the solution of the axial 
equilibrium equation and the determination of the fiber stress and shear 
stress. The inclusion of the matrix normal stresses in a modified failure 
criterion using the shear -lag model as well as using the coupled 
equilibrium equations of [6], with damage, is being considered by the first 


author and will be presented at a later date. 

The single equilibrium equation in the longitudinal direction is 


then: 

for all fibers, except N and N+l when y < L, 


W 

V 


d v 


n 


+ v 


dy 


n+l 


- 2v 


n 


v = 0 
n-1 


(18) 


for fiber N , y <_ L 


E 




d V,, 


G..t 2 . 
M ay 


” 1-1 


- - —x <y - = 0 , (19) 

g m 0 
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and for fiber N+l, y < L 


EpAyh 


A 


Jm. + v 

a 2 N+2 

dy 


v 4- "■*— x <v - &> ■«•• 0 

Vr ;G M . o/ 


( 20 ) 


Noting the coefficient of the second derivative term in the above 
equations, the following changes in the variables are suggested. Let 


e fV 

v CL,t 

v M 


dv 

n and a I = O a = E_ > then the normalized displacement 
t f 1 oo n F dy 

’ n 


y is defined by the equation 
n 


Ajh 

v = c, « V n • 


and the normalized shear stress r is given 


- 

t = n — t 
o 00 Epht o 


Algebraic manipulation then gives 


dV T E-ht dV 
i • . _ ' n _ o ? n 

fL " 00 <*n " f o Vp ’ 


(21) 


*l„ * ftsf (V n - Vl> * ^ (V n - W ' 


W 

V 


a , and l - 


EpV 1 

q-£— S , where 


M 


n, a , V , T , a and 8 are non-dimensional. 

’ n n’ o 

In these equations 2 ? , A ? , t, L and 1 are taken as actual fiber modulus, 
fiber cross-sectional area, lamina thickness and damage dimensions respectively 
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The quantities T q and G^/h are equivalent yield stress and lamina 
stiffness respectively and are to be determined experimentally. The yield 
stress, t q , should be reasonably close to the matrix yield stress obtained 
from a test using matrix material alone as long as the damage occurs in the 
matrix rather than along the interface or within the fiber. The quantity 
Gj^/b is felt to be less well defined as discussed above. 

The resulting non-dimensional equations are: 

For all fibers, except N and N*+l when n < ot, 


d 2 V 


an 


— — + y 

2 n+l 


2V + V . 
n n-1 


= 0 > 


( 22 ) 


for fiber N, n < a 


d 2 V, 


dn 


2“ ~ V N + V N -1 “ ^ <r1 - 3> = O'- , and 


(23) 


and for fiber N+l, n <_ a 


d 2 V, 


N+l 


dn 


V N+1 + V N+2 + ?<n - S > = 0 


(24) 


Defining a new unknown function f(n) such that 


f(n) = V N - v N+1 - f<n - 3> if n < a , and 


f(n) = o , 


n > a 
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with g(n) = V N - V N+1 for the same range of n values , the above three 
equations then become 


? 

d % 


n 


+ V - 2V + V ~ 0 


dn 


2 n+1 a n-1 


d 2 V, 


dn' 

d 2 V, 


-JL + v - 2V + V = -f(n)» and 
2 N+1 N N-1 ■ '■ 


(25) 


(26) 


N+1 


- 2V n + V M = f(n) . (27) 


dn 


2 N+2 N+1 N 


The^e differential -difference equations may be reduced to differential 
equations by introducing a- new function 


9) = =- Z V (n)cos(n0) from which (28) 

■' k o n 

n-0 , 


V^( n ) = “ /q V( n r 9 ) cos(n0 )d6 and ( 29 ) 


the three equations become 


2 - 

% f {—• - 2[1 t cos(0)]V}cos(n9)d9 

77 0 dn 2 


2 — 

— - 2 [l - cos(0)]V}cos(n0)d6 

77 0 dn 2 


,2W 


= 0 


(30) 


= -f(n)» fiber N, n < a, and (31) 


i- / 1T {^-2.- 2 [l - cos(0) ]V}cos(n0)d9 * f(n) , fiber N+1, n <_ a (32) 

77 0 dr ,2 


Making use of the orthogonality of the circular functions these three 
equations may be written as one equation, valid tor, all values of n 
and n» as follows: 
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ir^O “ 2 ^ 1 " cos( 9) }V}cos(n0)d9 = ~ <a-n> /* f(rO{cos[(N+l)0] 
dri . u 

- cos(N9) }cos(n0)d0 . ( 33 ) 

This equation is of the form 


2 tt 

-/ Q F(n,0) cos(n0)d0 V 0 for all n and n 

and noting the definition of V(n,0) in equations (28) and (29) it is seen 
that the function F(n,0) is even valued in 0 and therefore, if the 
integral is to vanish for all n, the function F(n,0) must be zero. 

The single equation specifying V(n,0) is then 


- 6 2 V = -<a - n>D 2 f(n), where ( 35 ) 

dn 

6 2 = 2[1 - cos ( 0 ) ] = b sin 2 ( 0/2) , and ( 36 ) 


D 2 = cos(N0) - cos[(N+l)9] . 


(37) 


It is very significant that the irregular boundary condition, ecuation 
( 1 ), of specified stress over a finite length, not coincident with either 
coordinate axis can be accounted for exactly and that the problem reduces 
to one differential equation which must satisfy boundary conditions alone 
the coordinate axes only. The ability to do so strongly depends on the 
form of the failure criterion. A condition in which both normal and shear 
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stresses were included generally would couple the axial and transverse 
equilibrium equations and yield a far more complicated set of differential 
equations. The apparent need to investigate such modifications is indicated 
by the results and, as mentioned above, is being considered. 

The solution to the problem of vanishing stresses and displacements 
at infinity and uniform compression on the ends of the broken fibers will 
now be sought. The complete solution is obtained by adding the results 
corresponding to uniform axial stress and no broken fibers to the 
following solution. 

The boundary conditions are then 

V Q * 0 as n 00 , ( 38) 

dV 

dn~ ” ^n = -1 ’ for 71 = °> broken fibers, and (39) 

V n = 0 , for n = 0, unbroken fibers. (ho) 

Using a technique such as variation of parameters to determine a 
particular solution to equation (35), the complete solution Satisfying 
vanishing stresses and displacements at infinity is 

V(n,e) = A(0 )e“^ n + J— sinh[6(n-t) ]f(t)dt (hi) 

where the unknown functions are A(9) and f(t). The remaining two 
boundary conditions give 
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dV (0) « _ 0 n 

— £ - f- {-6A(0) + D 2 /“ cosh(5t)f(t)dt}cos(n9)d0 = -1 (42) 

dri it 0 o 

for all broken fibers 

and 

2 

V (o) = - {A(0) - s inh ( 6 1 ) f ( t ) dt } cos ( n9 ) d0 = 0 (1+3) 

n it 0 o u 

for all unbroken fibers. 

Equation (43) is solved exactly by taking 

2. N 

A( 0) - t~ f* sinh(St)f(t )dt = 2 3 cos(m0) (44) 

6 0 m=0 m 

where the B are constants. Equation (42) then gives a system of N+l 

HI 

algebraic equations for the N+l constants in terms of f(n) which 

is, as yet, unknown. For the case of no damage the problem is then solved, 
i.e. see [7 ]. For example, consider the special case of no damage and 
one broken fiber. Then equation (42) gives 

— fl - 6 B cos(O)a 0 = -1 
it 0 o 

2B _ 

or / A 2 sin( 0/2 ) d0 = 1 , 

0 

Therefore 3 = v/8 = A(9) and 

0 



Ik 


V = fl A(0)e”^ ri cos(n0)d0 . 
n IT 0 


The maximum 


fiber stress is in the first unbroken fiber at n =0, and is 


* - weweue -1 
a dn it 0 3 

00 


or, for a 


fiber 


unit stress at infinity and an unloaded free end of the broken 


a 1 = h/3 ■ 


The normalized crack opening displacement, 2 V Q (0) , is ir/2. 

For matrix damage, a * 0, equation (42) must be supplemented by the 


condition that 


- f l s:n-e > , n < a , and 


g(ot) = T fl 


The constants 3^ and the function g(n) are then specified by 
requiring that equations ( U2 ) and (k5) be satisfied. Using equation (tl). 


dera'ole algebraic manipulation, the displacement of any 


and after consi 


fiber for all values of T| is 


I 
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. N 

v (n) = - f e^ n Z B c os ( m0 ) cos ( n0 ) d0 
n tt 0 nFO m 


+ f(t)( c n ( lt-n| ) - c n (t+n)}at , (46) 


where 


C n (5) = | /J f“ e " 6C cos(ne)ae 


Equation (42) then becomes 


N 

- fl (-6 E B cos 
tt 0 _ m 

m=0 


. D 2 e"^ t g(t)dt +• D 2 T a /g at} x 


x cos(n9)d0 = —1, n = 0, 1, ...» N 


(47) 


and equation (45) along with (46) gives 


. . N 

g(n) - — fZ e - ^ E B cos ( m0 ) { cos ( N 0 ) - cos[(N+l)6]}d0 
* ° m=0. . m 

+ |/o g(t){c N (|t - n|) - c N (t + n) - c N+1 (|t - n|) + c N+ i (t + n)}dt 


"2 ^3 - n| 


The condition that 


- c, T ( t + n) - c N+1 (|t - nj ) + c N+1 (t.+ n))dt . 


N+l 


>8) 


g \ C L j — T 0 


(49) 
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also must be satisfied . 

Physically, it would be more direct to specify the applied stress, 
and the number of broken fibers, N, and determine the damage zone a and 6 
depending on given yielding and splitting conditions. As a and 3 appear 
in the limits of the above integrals this is not convenient mathematically 
and it is easier to specify the number of broken fibers, N, and the damage 
zone a and 3, and compute the required applied stress o ra . 

These equations were solved as follows : 

I. An initial set of constants B was determined for the problem of no 
damage, ct = 3 = 0 in equation (47)> i*e., 

N * * . ,, , 

2 B 6 cos(m8)cos(n0)d0 - 1 , (50) 

n m ir 0 
m=0 

n = 0 , 1, . . . , N 

II. These initial constants were then substituted into the integral 
equation (MB) and, along with equation (4$), the function g(n) 
and were determined using the desired values for ot and 3. 

lit. Using g(n) and f dJ a new set of constants, 3 m , was computed from 
equation (47) with the desired values cf d and £< 

XV. This procedure was repeated until the unknowns changed less than a 
prescribed amount with additional iterations . 

In the above solution the Unknown function, g(n ) , was assumed to 
be piece— wise linear ever the interval 0 < n . <. ot of the form 

g*(u) » y't + Ytu, i = 1,2,. . .k 

V 
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when the interval was divided into k equal subdivisions. The function 
g(r|) then contained 2k unknowns with one additional unknown being T 0 . As 
g(q) is the displacement difference it should be a positive, monotonically 
decreasing function and its representation as a piece-wise linear function 
should be sufficiently accurate. The (,2k + I i equations were obtained by 
requiring that the integral equation, equation lk8 ) , be satisfied at the 
(k + l) end points; (k - l) equations resulted from the requirement of 
continuity of the function g(n ) between adjacent intervals and the last 
equation was given by g(a) => T q . 

With the longitudinal displacement v^ now known the transverse dis- 
placement u q is obtained by solving equations (4), (6), and (8). 

Equation (10 ) gives the matrix normal stress in terms of u^. This solution 
is recorded below for completeness. 



+ sin [(N 


+ 



sin ( n6 ) d8 
1 - cos (0) 


(51) 



SOLUTION AND RESULTS 

The numerical solution to equations (Vf), (k8) and ( U9 ) was developed 
as described above and was solved using the NASA-Langley Research Center 
CDC computer. Convergence was fairly rapid with the number of subdivisions, 
k, required in the representation of the function g(n) being on the order 
of fifteen and the number of iterations necessary to give a change in the 
normalized shear, f, such that (? i+1 - < ICT^, on the order of ten. 

Larger values of the damage region* ct and 3 , required an increase 
in both the n umb er of subdivisions and number of iterations, with a, S < 2 
being relatively small values and ct, $ > 6 requiring increased accuracy. 
Computation time varied bath with, the above parameters and with the 
number of broken fibers with typical times being on the order of one minute 
fpr ct, 3 ~ 2 and N = 1, to twenty minutes for a, 3 « 4 and N = 15 • 

One problem of particular significance is the behavior of a laminate 
after damage develops and the investigation of the potential for continued 
longitudinal yielding or splitting or for notch extension due to progressive 


fiber breaks. Some typical results are given in Figure 3 for one and seven 
broken fibers using a two /one split strain— to— yield strain condition as show n . 
The maximum fiber stress normalized by the laminate constant 
is plotted against the normalized applied stress where the maximum 
fiber stress is always found to occur in the first unbroken fiber at the 
end of the split, (n =3). 


These results are representative of many different cases worked, i.e. 

increasing the strain ax which splitting occurs simply increases both O ^ 

and a at the ooint of splitting but the nature of the behavior is 
00 “ 


unchanged. 'That is, in all cases once the split forms the fiber unloads 
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and the split length becomes unbounded under a five to ten percent further 
increase in applied stress, at which time the behavior reduces to that of 
an unnotched laminate with * a . The net section fracture stress is 
then independent of initial crack length as [ 3 ] • 

This predicted behavior indeed happens it. unidirectional graphite/ 
epoxy or boron/epoxy composites. However, for a less brittle matrix 
material such as aluminum some small splitting has been observed [3] and 
[7], but it is stable and the stress in the fiber continues to increase 
with increased remote stress until the ultimate fiber stress is reached. 

The model then appears to give reasonable results for large splitting but 
does not account for small, stable longitudinal splitting. It is of 
considerable interest to determine if the inability of the model to 
predict the behavior of stable splitting is due to the assumed failure 
criterion or the shear— lag model. As mentioned above, a more complex 
failure criterion including the matrix normal stress as determined from the 
shear -lag solution as well as more complete shear stress-displacement 
assumptions are being studied. 

Surprisingly, in view of the above difficulties, the present model 
approximates the behavior of boron /aluminium amazingly well for strains 
such that splitting does not occur, and gives an accurate estimate of 
the laminate strength as a function of number of broken fibers (crack 
length) as well as crack opening displacement. The following results 
are for no splitting (3=0) and conclusions for a specific boron /aluminum 
laminate will be drawn from the general results. 
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Figure !*• gives a plot of the maximum fiber stress as a function of the 
applied stress for different numbers of broken fibers . Defining a function 
S equal to the normalized ultimate fiber stress as, 

s = gult 1 °^ 

s '= vV- ’ 

the required applied stress for a specified ultimate fiber stress is 
then given by the intersection of the horizontal lines corresponding to 
specific values of S as shown. 

As the applied stress is increased the matrix undergoes plastic 
deformation as indicated by Figure 5. For all ranges of a indicated, the 
maximum shear strain which is at ry = 0 between the N and N+l fibers , 
is no greater than three times the yield strain and therefore for a boron/ 
aluminum composite splitting probably would not occur. The lines of 
constant 3 values on Figure 5 are obtained by locating the corresponding 
(o^jN) points from Figure U. Figure 6 is also developed from Figure 4 
and gives the strength of the laminate as a function, of number of broken 
fibers (or crack length). That is, a is the remote stress required to 
give a particular ultimate stress in the first unbroken fiber. 

In addition to these results the displacement of the center broken 
fiber, or therefore one-half the crack opening displacement, is obtained 
and is depicted in Figure T as a function of applied stress for different 
numbers of broken fibers. As the number of broken fibers (crack length) 
increases it is seen that the matrix yeilding contributes a constant 
proportion of the total crack opening displacement. This is consistent with 
Figure 5 where, if the damage length, a, is plotted against number of broken 
fibers for constant remote stress , “Figure 3 is obtained and the relationship 


is linear as shown. 
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Using these general results the values of x • and the stiffness constant 

G /h will now "be determined by comparing, with the experimental study of 
M 

Averbuch in 18] for unidirectional boron/ aluminum. The laminate used in [8] 
had the following material and geometric properties : 

Ep = 475 x 10 9 Pa, 

Ay = 1.59 x 10 -8 m 2 , (D * 0.1427 mm), 

t = 0.165 mm/ply, eight plies, 

a , . = 3.98 x 10 9 Pa, 
ult 

w = width = 25.4 mm, and 
fiber centerline spacing = 0.178 mm 

For a laminate having seven broken fibers, which corresponds to a crack 
length of about 1.27 mm, the load vs. COD curve is given in Figure 9. An 
approximate "best fit" of this curve with [8] requires that the yield stress 
and stiffness be 

Q 

x = O.IO 9 x 10' Pa, and 
o 

G M /h * 65.4 x 10 12 N/m 3 , 

from which the normalized ultimate stress, S, is 


= |Vf , i, 1 

t o yv* 

For reference, the damage length is then L - 4.71 ot 

laminate constant T = 3.65 T . 

o o 

'lew , using these values, the corresponding curve 
a crack length of 5 mm is plotted in Figure 9 and the 


fiber spacings and the 

for 29 broken fibers , or 
comparison is seen to ce 


very good. 
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Further, referring to Figure 6, the failure stress of the laminate as measured 
in [8] is seen to be of the:same form as the predicted value, although the 
model gives a much larger decrease in strength for small crack length than the 
experimental result. The value of S determined from the COD measurements 
does, however, predict the failure stress reasonably well for longer crack 
lengths, The significant point is that the simple shear-lag model does relate 
crack opening di s plac ement and laminate strength. Changes in the shear-lag 
assumption and/or the failure criterion may well improve the agreement with 
experimental results , and an investigation into which modifications are 
important should lead to a better understanding of the fracture process . 

It should be noted that the value of T q given above is close to the 
measured value for a homogeneous aluminum specimen and is also approximately 
equal to the value determined experimentally by Peters [3] for a boron/ 
aluminum compos ite . The equivalent stiffness, G.^/h, for h — 1, 18 x 10 m, 
which is the center— line distance between fibers , gives a shear modulus 
G = H.6 x 10.9 Pa. This la on., the order of one half the value for aluminum. 
Smaller values of h, of course, give smaller values of G^. 
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CONCLUSIONS 

It has "been shown that a model based on the shear-lag assumptions and 
using a shear stress failure criterion of an elastic-perfectly plastic 
matrix material is capable of predicting the behavior of a unidirectional 
boron /aluminum composite quite accurately foi strain levels below those 
sufficient to cause splitting. For brittle matrix composites the model 
gives equally good results for both yielding and large splitting. 

The fracture stress is then predicted for the two limit cases of large 
yielding with no splitting and large splitting. For yielding alone the 
fibers nearest the notch continue to increase in stress as the applied load 
is increased and the fracture sbresis is crack length dependent. In a 
brittle matrix longitudinal cracks develop tinder low loads and become 
unbounded for a small increase in load, thereby negating the effects of the 
notch and the fracture stress is crack length independent . 
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Figure 3. Yield length as a function : c a -=lied stress, 
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Figure 6. Applied stress as a function of number of broken fibers for 
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Figure 8. Damage length as a function of number of 
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